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Abstract 

The problem of coupled Fermi-Bose mixtures of an ultracold gas near a narrow Feshbach res- 
onance is approached through the time-dependent and complex Ginzburg-Landau (TDGL) theory. 
The dynamical system is constructed using Ginzburg-Landau-Abrikosov-Gor'kov (GLAG) path in- 
tegral methods with the single mode approximation for the composite Bosons. The eguilibrium 
states are obtained in the BEC regime for adiabatic variations of the Feshbach detuning along 
the stationary solutions of the dynamical system. Investigations into the rich superfluid dynamics 
of this system in the shallow BEC regime yield the onset of multiple interference patterns in the 
dynamics as the system is quenched from the deep-BEC regime. This results in a partial collapse 
and revival of the coherent matter wave field of the BEC, whose temporal profile is reported. 

1 Introduction 

Over the course of the last two decades, there has been a groundswell of theoretical and experi- 
mental interest in ultracold gases of alkali metals confined in optical and/or magnetic traps. These 
systems of ultracold gases have proved to be extremely robust and tunable systems for studying 
condensed matter physics in regimes that are inaccessible in solid state systems ll|]. During this 
time, physicists also became specifically interested in condensates of Fermionic alkali atoms (such 
as 6 Li), and obtaining a BCS superfluid of Cooper pairs similar to those seen in solid state supercon- 
ductors (BCS theory). If the effective attraction between Cooper pairs can be rendered sufficiently 
strong, Cooper pairs of Fermions are no longer merely correlated and far apart as in traditional BCS 
systems, but have much smaller correlation lengths approaching the interparticle spacing. Thus, 
they can be treated as composite Bosons, causing the system to undergo a BCS-BEC crossover to a 
BEC superfluid. 

Theoretical work on the single channel model of the BCS-BEC crossover over the course of 
the late 1980s and 1990s [|2H_6] motivated numerous experiments with laser cooling and trapping 
of Fermions, culminating in the observation of this crossover in the early 2000s |7- 10]. In the late 



1990s, the more informative dual-channel model of resonant superfluidity was proposed, by Holland 
et al II 1 lfl . and by Timmermans et al II 1 211 - This model, built from the Timmermans' Hamiltonian, 



is a generalization of the Dicke (Tavis - Cummings) model in guantum optics 111 311 - The attractive 
interaction between the Fermions can be controlled by tuning a homogeneous magnetic field to their 
Feshbach resonances [12], which are caused in a two-particle system by coupling bound states in 
a close channel with states in the scattering continuum. The effective scattering length can be 
tuned simply by varying the external magnetic field that controls the net magnetic moments of 
the different channels. In this dual channel model, sufficiently strong resonances cause bound 
states to form in the closed channel. Thus, Cooper pairs of Fermions can physically combine into 



Bosonic molecules Ill4f l, entities referred to variously as composite Bosons, diatomic molecules, or 



quasimolecular states in the literature. These composite Bosons have extremely long lifetimes near 



a Feshbach resonance [15, 16], and repel each other |17], facilitating their condensation to a BEC 



superfluid H 1 8 . This area of research has gained enormous interest over the last two decades. 
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An overview of the phase portrait of a population balanced Fermi gas can be found in dob, and 
reviews of the current status of research can be found in J21H22I1 . 

The possibility of observing the nonequilibrium dynamics of coherent quantum states via quench- 
ing (diabatic variations of the system parameters) is one that is unique to systems of ultracold 
gases, and unavailable in other similar condensed matter systems. Post-quenched dynamics in cold 
atom systems have been studied and reported. In particular, the collapse and revival of a post- 
quenched coherent BEC state in an optical lattice 112 3H generated considerable interest. Quenched 
dynamics have also been investigated for Fermi-Bose mixtures, experimentally 12411 , as well as theo- 
retically lfl3l.[26il30ll. The possibility of observing a collapse and revival phenomenon as the system 
is quenched past the BCS-BEC crossover has been raised Il32n . More recently, nonequilibrium dy- 
namics in condensed matter physics have also been described by the complex Ginzburg-Landau 
equation [33]. The nonlinear damping in the dynamics produces very rich and diverse behavior 
in Fermi-Bose mixtures I134n . The mean field dynamics of BCS superconductivity have been ob- 
tained from microscopic models via Ginzburg-Landau-Abrikosov-Gor'kov (GLAG) theory. A similar 
approach has proven successful for BCS-BEC systems as well, thus leading to their description by 
the TDGL equation for the single channel case [4]. More recently, the applicability of TDGL dynam- 
ics have been demonstrated in the dual channel case I134n . This motivates the use of this treatment 
to study the dynamics of coherent matter waves in BCS-BEC systems in this report. 

This paper focuses on the dynamics of the TDGL equation in coupled BCS-BEC systems, and 
the dynamics of quantum quenching therein. The most general state of this system across the 
phase diagram is that of a Fermi-Bose mixture, where the composite Bosons coexist with correlated 
Fermions, and their dynamics are linked in the mean field by the two-channel scattering process 
through coupled Ginzburg-Landau and Gross-Pitaevski equations. The mixture is characterized by 
two distinct phases, the Fermi (BCS) superfluid phase given by the order parameter and the 
Bose (BEC) superfluid phase given by the order parameter The Fermi superfluid consists of dis- 
tinct Cooper pairs, and the Bose superfluid consists of quasimolecular Bosons. This 2-channel model 
is more faithful to the microscopic nature of the system than the single-channel model, especially 
when Feshbach resonances are involved. The dynamics of the coupled BCS-BEC phases is richer 
than that described by a single-phase model in the single channel case, since the Ginzburg-Landau 
dynamics of the latter is only the damped one except in the BEC regime. This is especially true if the 
Feshbach resonance is narrow, which is the case being discussed in this work. Though many current 
experiments correspond to broad resonances, the consideration of a narrow Feshbach resonance 
here is not unrealistic at all |35, 3^]. The paper reconstructs the time dependent Ginzburg Lan- 
dau equations for the dynamics of a Fermi-Bose mixture in the two channel case starting from the 
many-body functional field integral for the Timmermans' Hamiltonian. This Hamiltonian describes 
the Fermions by the BCS Hamiltonian and the composite Bosons by the Bose-Hubbard Hamilto- 
nian |12|], and includes a resonant coupling between the two species. The path integral is written 
as a functional of the superfluid gap parameter A and the boson coherent state (all of them are 
taken to lie in the zero momentum state) bo, taken to be a c-number. The mean-field dynamics is 
subsequently derived, and the equilibrium densities evaluated at the stationary solutions. The pa- 
per then details investigations of the dynamics of quantum quenching, and reports the possibility 
of collapse and revival in the full matter wave of the BEC at large times. Section [2] begins by outly- 
ing the formalism that obtains the TDGL dynamics from the Timmermans' Hamiltonian. Section [3] 
reports the study of the stationary solutions of the TDGL dynamical system in 3 dimensions. The 
chemical potential and condensate fractions are evaluated as the Feshbach detuning is varied adi- 
abatically. Section|4]looks at the dynamical evolution of the BEC (which can be seen directly in the 
lab via time-of-flight absorption) as the Feshbach detuning is varied diabatically (ie 'quenched') to 
the shallow BEC regime. Concluding remarks are made in Section[5j 



2 Dynamical Equations of Motion 

The treatment that obtains the coupled TDGL dynamics of this system closely follows that which 
obtains the conventional TDGL dynamics of the single channel model, as applied by Huang, Yu and 
Yin 13711 - The formalism is applied to the dual channel case in a manner similar to the treatments by 
Machida and Koyama [34], except for the nature of the Bosonic states, which is approximated by a 
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single mode. It is assumed that, at T = 0, the time-scales of the dynamics are sufficiently weak so 
as to not induce Boson excitations above the ground state, an assumption justified in the context of 
dynamics in greater detail in the literature 1 26]. This approximation greatly simplifies the dynamics 
by removing any spatial information in the composite Boson amplitudes from the beginning. This 
approximation also allows for rapid transitions near unitarity. 

The D— dimensional Fermion and zero momentum composite Boson fields in this system are 
represented by the operators (j) a (x) and 60 respectively, with the index a representing the 

Fermion pseudospin. The dual-channel Timmermans' Hamiltonian 111 211 H tm for a Fermi-Bose mix- 
ture at T = for a unit volume is 

H tm = J A D x x H tm (x), (1) 

where 



H tm (x) = \ #U X ) l h ( r ) ~ Mf] &r(ar) J - |u J ?|^|(a;)^|(x)^(a;)^ 1 .(x) 

+ [2v - fi B ] bl(t')b (t') + u B blb (blb - l f 

+ 9r bl{t')(t>^{x)(j)^(x) + h.c 



(2) 



Here, x = (r,f). Eguation [2] describes a system of ultracold electrically neutral two-component 
Fermions interacting attractively. The first line in eguation [2] represents the Fermi-BCS part of the 



Hamiltonian. Here, h(r) = — j— + V(r) 



is the single particle Hamiltonian, and m is the Fermion 



mass (the mass of the composite Bosons is thus 2m). The second line represents the Hamiltonian of 



the composite Bosons 112 



32h . Here, the Feshbach threshold energy (also called the Feshbach 'de- 



tuning' from the molecular channel to the continuum B12I1 ) is represented by 2v, and ub represents 
the amplitude of the repulsion between the composite Bosons. The final line describes the Feshbach 
resonance that leads to the Fermions binding to (or dissociating from) the composite Bosons, with 
g r representing the atom-molecule coupling. Note that the chemical potentials satisfy (j,b = Z^f- 
The path integral grand partition function Z is defined by 



Z = 



]D[&*,&]e _s *V 



(3) 



where D[<f>, cf>] and D [&*,£>] are the path integral measures of the Fermion and Boson fields respec- 
tively, and t' is the imaginary time. The action S^b is given by 



'#0 



Mx)dt>Mx) + b* (t')dt>b (t') + n tm $A,b* ,b ) 



(4) 



This integral can be evaluated by introducing the macroscopic gap parameter A(t') into eguation[3] 
via the Gaussian identity 

.A*(i')A(i') 



J V[A*, A]exp - J At 



= 1, 



(5) 



where the Fermion field (characterized by the gap parameter) is spatially homogeneous due to it's 
coupling to the spatially homogeneous Bose field b by momentum conservation. Performing the 
Hubbard-Stratonovich transformation, A — > A + \up |<Al.</>t an d A* — > A* + cancels out the 

four- Fermion term from H tm in eguation |4l Now, integrating out the Fermion fields 4>, <j> remaining 
in eguation [3] using formal Grassman calculus yields 



D[6*,6]D[A*, A]e 



-SAb 



(6) 



where 



S A b - / dt'l [2u - fi B } \b (t')\ 2 + u B \b (t r )\ 2 [|6 (<')| 2 - l] + b*(t')d t ,b (t') + 



\Mt')\ 2 

\u F \ 



J d D xx At' lndetM(x). 
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Here, 



M(x) 



d t > - tip + h(r) A(t') + g r b Q (t r ) 
A*{t')+g r b* (t') d t ,+ii F -h{r) 



The action Sa& is split into a field independent part So = J d D xdt' lndet M (x), where 



Mo (a:) = 



d t > - fi F + h(r) 

d v + I-Lf - h(r) 



(7) 



(8) 



and a field dependent part 5 e // which vanishes when A and b do. Expanding S e ff to the fourth 
order in A + g r b 1 37] and performing gradient expansion results in 



S eJf « / d v x\ d [A*(f ) + ft' [A(i') + 3r 6 (t')] + 



|A(t')| s 
k/ F | 



Here, 



o|A(0 + g r b Q (t')\ 2 + ^b\A(t') + g r b (t')\ i + 

[2u-^ B ]\b Q (t')\ 2 + u B \b (t')\ 2 [\b (t')\ 2 - 1] +b* (t')d t ,b (t')\- O) 



& v x' x Q(x-x'/2,x + x'/2), 



b = 



D 

/ JJd^Xj x R(x,xi,x 2 ,x 3 ), 



the coefficient d is obtained from [34], 



d = Um / d D x' x 



Q(x-x'/2,x + x'/2). 



(10) 



(11) 



Here, d^x = d-°x x eft', where d^a; is the measure of integration over all D spatial degrees of 
freedom, and 

Q(xi,x 2 ) = G + (xi,x 2 )G-(x 2 ,x 1 ) 
R(xi, ...,x 4 ) = G + (xi,x 2 )G-(x 2 ,x 3 )G + (x3,x 4l )G-(x4,,xi), 



with x' = (r',t'), and the Gor'kov Green's function, G(x) = M (x) = 
obtained from the Green's function of a noninteracting Fermi gas i.e. 

[dv TV± h(r)} G±(x - x x ) = 6 D (r - r x )5(t' - t' x ) 



G+(x) 
G-(x) 



(12) 
which is 

(13) 



Note that the product in the expression for b above is over the spatial dimensions D only, whereas 
the measure is over the space-time dimensions V. The mean field eguations of motion of the order 
parameters can be obtained by eguating the functional derivatives SS e ff/SA*(t') and SSeff/Sb^t) 
to after analytically continuing the time to the real axis by substituting it —> t. This yields the final 
dynamical eguations for this system 



*i + 17 (*i - * 2 ) - ia&\ + z/3 1 \1> x | 2 \I> ! = 
* 2 + 2i\V 2 + 2ix|*2| 2 *2 - i«7(*i - * 2 ) = 0. 



(14) 



Here, 



*2 



A + g r b Q 



\»f\Vn 

9rb 



(15) 
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where J\f is the total (Fermion) particle number. The constants expressed as Greek letters in equa- 
tion [14] are given by 



a\/j. F \, A 



IMF I 



\fx F \d, 



' |P 1 ^ ' (16) 

7=|£f|, X = du B \fx F \a 2 Af. 

k = g 2 d, 

Finally, note that time has been rendered dimensionless via the transformation t — > n^n x t> an d 
the chemical potential is presumed to be negative. Thus, the dynamics of the Fermi-Bose mixture 
is that of a system where the Fermi and Bose fields evolve according to coupled Ginzburg-Landau 
Gross-Pitaevski-Bogoliubov dynamics. The coupling is caused in the Ginzburg-Landau case by the 
order parameter A getting nonlinearly dressed by the Bose field g r b , and in the Gross-Pitaevski 
case by a harmonic coupling to A. 

It is also noted that the dynamical system bears a resemblance to that of a molecular BEC of 
atomic Bosons in a resonance effective field theory as proposed by Kokkelmans and Holland in 



2002 [25]. The role of the pairing field of noncondensed atoms (represented by the 'anomalous den- 
sity' of noncondensed pairs) in that system is assumed by the Fermion gap parameter A (related to 
the anomalous Cooper pair density by A* ~ 2~Zk( ak t a -kj.)) in this one. This strengthens the analogy 
between Cooper pairs in Fermi-Bose systems and noncondensed atoms in Bose-Bose systems. 

The dynamics is investigated for a system confined in a 3 dimensional box where the confinement 
V = C 3 is much larger than the inter-particle spacing, effectively treating the trap as homogeneous. 
The constants a, b and d can be evaluated from equations [TO] and [TT] by using the Green's function 



for a free particle. This yields 134. 13711 



a — Z)|k|<fc R 2e fc ' & — S|k|<fe fl 4i|' d ~ J2\k\<k R 4iJ" (l 7 ) 

Here, = ^ — \i F is the Fermion energy and k R ■ — ^| is the renormalization cutoff in momentum, 
placed to counter ultraviolet divergences. An important caveat here is that the integrals in equa- 
tion [T7] contain singularities if [i F > Thus, this formalism breaks down in that regime, which 
corresponds to regions where the BCS state dominates over the BEC state, well beyond the current 
region of interest at a negative fj, F . Going to the continuum limit by substituting for the formal sum 
J2 J d 3 fc in equations [T7] and performing the integrals results in 

mkjiEp (2me F ) kn 
a[e F ) = - - arctan ■ 



2tt 2 4tt 2 V^eF' 



qi \ (2me F ) 3/2 
P{eF) = ^28^ XA ^' 

e F mk R e F e F 

. gl (2me F f 2 

K ( £ F) = ^T-X ~2 . 

K (cf) = x 



(2me F ) 3/2 



er(e F ) 



327r e F 
£F 



9r 

1 \ u B a 2 (e F ) (2me F ) 3/2 

x(£f) = x x Af , (18) 

327T e F 

where the dimensionless constants in Greek letters are now functions of the chemical potential 
e F = I^fQ- and the ^-dependence on A„ has been emphasized by subscript. In the equations above, 



1 Note that, in general, e F as it is denned here does not equal to the Fermi energy. 
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the relation 

1 _ J_ X ^ J_ 

1 " 1 Fl |k|<fc R k 



(19) 



has been used to obtain the expression for 7. Here, e° = u° F is the bare interaction 47ra s /mV, 
and a s is the s-wave scattering length controlled by Feshbach resonance. This relation comes about 
as a consequence of renormalizing the BCS gap equation so as to counter ultraviolet divergences 
in the single channel model [37|] . Finally, note that the inability to renormalize this theory 

analytically disallows taking the theoretical limit of kji — > 00 for some of the constants above. This 
problem will be discussed further in the subsequent section. 



3 Stationary solutions and Chemical Potential 

For a complete phenomenological description of the dynamics, number conservation has to be sat- 
isfied and used as a constraint at t = to obtain the chemical potential tp = \fip\ where (ip = ^ . 
In order to do this, the action S e ff from equation[9]is investigated at small temperature j3 = l/k B T. 
Since we only consider regimes with negative chemical potential, all the Fermions in the gas are 
correlated and no 'free Fermions' remain H. Thus, in the stationary case, 



SefM 



at' 



Up 



a|A + g r b Q \ 2 + [2v - fi B ] \b Q \ 2 + u B \b \ 2 



1 



(20) 



where the quartic contribution has been neglected. The temperature dependence of all constants 
can also be neglected and the expression above simplified to get 



u F 



- o|A + g r b \ 2 + [2v - mb] |6 | 2 + u B \b \ 



The Helmholtz free energy at T = is calculated from 

n = lim — -lnZ(/3), 

P— s-oo p 

where 

Z{fi) = J D[b*,b}D[A*,A}e- S 'ff^l 
Simplifying equation [23] by taking the mean field (ignoring fluctuations) yields 

Thus, the Helmholtz free energy Q at T — is given by 

IAI 2 



Up 



a\A + g r b \ 2 + [2v - N + usN 



Imposing number conservation by using Af = —dfl/dfip, 

M=^\A + g r b \ 2 + 2\bo\ 2 . 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



where fj, = fip = —ep, and the dependence of a on fi is obtained from equations 1171 Writing the 
equation above in terms of dimensionless variables, 



where 



2<r 2 \<f 2 \ 2 +£1*1 1 
, da 



1. 



e = 



9/i 



(27) 



(28) 



2 The number density of uncorrelated Fermions, obtained from the 'free Fermion' term So = J d D xdt' lndet Mq(x) in 
section[2] is = @(^f — e^), which is if [if < 0. 
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Equation [27] is the constraint that fixes the number of particles M via the chemical potential 
Furthermore, in the stationary case, equations [14] reduce to 



7(*i 



*9 



■/3|*i| 

h - 



9i = 

2) = 



0, 
0, 



(29) 



where $2 are the stationary solutions of the dynamical system. The trivial stationary solutions 
of the dynamics in the ^1, ^2 phase space, ^1,^2 = 0, are realized at temperatures above the 
critical temperature T c 113 711 . The nontrivial ones are the locus of points satisfyinq the equations 
above plus the chemical potential eguation l27l Thus, the 3 unknowns ^1, ^2, £f are solved from the 
3 simultaneous equations 



7 (e F ) (#1 - tf 2 ) - a (e F ) tf x + j3 (e F ) |* x | 2 # x 



2A„ (e F ) * 2 + 2x (e F ) |tf 2 | 2 *2 - K (e F ) 7 (e F ) (*i 



* 2 = 



2a 2 (e F )|* 2 | 2 +£ 2 (e F )|*i| 2 -l 



(30) 



As explained in section [2] the BCS-dominant reqime is not entirely accessible in this formalism. 
However, the predominantly BEC qapless reqime can be obtained from eguations l30l In that reqime, 
A = ie ^1 = #2- I n the case of noninteractinq Bosons, ie ub — (x = 0), the second equation 
of [30] necessitates that A„ = ie e F = — v, in aqreement with mean field results H32n . This is also 
consistent with the physics of the system, since the chemical potential is the enerqy required to 
remove one Fermion from the system. In the qapless and noninteractinq BEC dominant case, the 
majority of the Fermions are dimerised with bindinq enerqy — v , and a dimerised Fermion requires 
v enerqy to dissociate from the dimer and free itself. For the more qeneral case, 



x 



a 



Clearly, \ v and so v has to be neqative for this to be true (\ is always positive). Also 



Q 

1 



K 

x 



0. 



(31) 



(32) 



Since fc F 3> I, so is a. Therefore, for the equation above to hold, X u and therefore v must be larqe 
and neqative for a qapless BEC dominant state to exist. This criterion is in aqreement with the 
physics of the system. For larqe neqative v, the bindinq enerqy and therefore the molecular affinity 
of the Boson dimers will be larqe and neqative, facilitatinq the dimerisation of the majority of the 
Fermions [14]. In order to calculate the chemical potential in this reqime, the relation a = f3\^/ 2 \ 2 is 
substituted in equation [32] Then, equations [18] are applied, yieldinq 



e F 



9r\b \ 



u B e F 

Solvinq this quadratic equation in e F and rejectinq the unphysical root, 



1 



e F 



1 



1 + 



1 



r UB\bo\ 



1 V 
- 2 U B) 



1/2 



V-~U B 



g r u B \bo\ 



(33) 



(34) 



In the case of noninteractinq Bosons, equation [34] qives e F = — v as above. Furthermore, in this 
reqime | fc» 1 2 = AT/2 (no BCS state, all Fermions are dimerised), or |^i.2| 2 = l/2er 2 (e F ). This sim- 
plifies equation [27] to £ = or da /dp = 0. Fiqure [T] shows plots of da/ dp as functions of v (where 
v = pf) for m = 1, V = 1 and several larqe values of fc F . Note that da/dp and therefore £ indeed 
vanishes in the limit v — > —00, which is where the BEC reqime is expected. 

Continuinq with the case of noninteractinq Bosons {ub = ie \ = 0), equations [30] reduce to 

7 (e F ) (t x - * 2 ) -a(e F )*i+/3(e F )|*i| 2 *i = 0, 
K (e F )*2 - K(e F ) 7 (e F ) (§ x - * 2 ) = 0, 

2a 2 (e F )|* 2 | 2 +e 2 (e F )|*i| 2 -l = 0. (35) 
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with X v (e F ) 



3/2 



327T 



[1 + ^-]. Solving the first two equations of l35l vields 



*2 



- (a - -yrj v ) 
(l-»7v)*i, 



1/2 



JO 



(36) 



where r\ v 



The problem of renormalizability of this dual-channel theory mentioned in 



i+(a„/k 7 ) • 

the previous section manifests here. The inability to renormalize happens because of the quantity 
a — 7J/V in the equations above. In the single channel case, when g r ,n —> and t] v — > 1, this 
expression converges as k r — > oo, as can be seen using equations [18] yielding 



lim 



lim (a — 777) 



(37) 



However, finite nonzero g r leads to rj vanishing as &r — > 00, yielding limfc R _>. 00 (a — 777) ~ a which 
diverges for arbitrarily large kn (see equations 1 18l l^l. 

The difficulty described above can be worked around by noting that the allowed momenta need 
to have a finite cutoff k^, which is governed by the physics of the system. For instance, in the case 
of solid state BCS systems, the cutoff is given by the Debye frequency on account of the physical 
origins of the electron-electron attraction in the slowly relaxing lattice vibrations. Ordinarily, the 
choice of cutoff is decided by the range of the interatomic potential, with k^ ~ S^v/gf if the 
bare repulsion between the Fermions is neglected [12, 38]. However, for small confinements and 
small detuning, the maximum momentum allowed in these cold atom systems cannot be greater 
than 2-k/R, where R is the size of the atoms. Experimentally, these atoms are very cold, and never 
move so fast as to cause inelastic scattering beyond the Feshbach resonance, since such excitations 
will dynamically alter the internal degrees of freedom of the atoms themselves. Thus, this choice 
of ku keeps a reasonable upper bound in all of the integrations. For Fermionic atoms like Li, R 



is 0(10 



-io„ 



which needs to be rescaled with respect to the unit length viz. the confinement 



dimensions. The confinement, however, is highly tunable in cold atom systems, with a fairly wide 
range of permitted values. Thus, using a wide variety of k^'s to get quantitatively correct results 
should be permissible, once the trap sizes are adjusted accordingly. In this report, Hr are chosen to 
be O(10) in units of inverse trap size. These should be attainable in small volume systems of O(10 2 ) 
atoms (obtained experimentally via laser culling by Chuu et al 1 39]) confined to length scales of 
O(10 _9 m) (obtained experimentally by confining ultracold gases using atom chips [40]). 

Equations [36] are nothing more than the gap equation for the BCS-BEC system in position space. 
In all these equations, a, f3, 7, n, X v , rj v are functions of e^- Note from the above that if k — > ie 
g r — > 0, then r\ v — > 1, and ^ 2 — > 0, which will give a pure BCS until ep goes to a regime where a can 
no longer be evaluated without running into singularities (as per section 0, causing this formalism 
to break down. Plugging the values of $12 above to the final equation of l35l vields 



1 



1 



a(e F ) - "t{e F )r] v (e F ) 



|2a 2 ( ei ,)[l - Vu (e F )] 2 -K 2 (ef) j- (38) 



Equation [38] is a transcendental equation and needs to be solved numerically for jj, F (—e F ). Fig- 



3 It should be noted, however, that this divergence goes away at T ^ due to the use of thermal Green's functions 
(instead of the ones from eguation ll3> in the evaluation of the constants in equations llOl |34[l l37ll . 
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Figure 1: (Color online) Plots of da/dfi as a function of \fj, F \ for kR — 100 (red), 500 (green), and 10 9 
(blue). Here, m = 1 and V = 1. Notice the rapid rate of convergence as fc# — > oo. 

ure [2] contains the numerical results of evaluating the condensate fractions and chemical potentials. 
These have been obtained by setting Af to 100 and using representative values of the interaction 
parameters. The chemical potential fip was obtained by numerically solving equation 1381 using 
Newton-Raphson methods. The default working precision was kept to the eighth place of decimal. 
Figures [2] (a), (b) and (c) contain plots for g r — 25, and figures [2] (d), (e) and (f) contain plots for 
g r = 40. Figures [2ta) and (d) are plots for (i F /i> as a function of v. Note from these figures that, 
for sufficiently large —v, [i F ~ v in accordance with the analytical results in the BEC regime. Fig- 
ures 0b) and (e) are plots for the condensate fractions n F (Fermions in BCS) and (molecular 
composite Bosons in BEC). The fractions are computed after solving equation 1381 and substituting 
the corresponding values of | ^ 1,2 1 2 (obtained from equations [36) into the relations (see equation !27H 



np 

»B 



e\*i\ 2 , 

a 2 \^\ 2 . 



(39) 



The rise in the BCS superfiuid density can be clearly seen in figure [2b) around v ss —50, and 
in figure Oe) of v ss —150. These results agree qualitatively with more sophisticated theories of 



BCS-BEC systems 114111 . considering that numbers are highly sensitive to the choice of the renor- 
malization cutoff kR. In the case of interacting bosons, equations [30] are simplified by operating in 
a regime where the chemical potential is weak enough so that (3 can be ignored. This yields 



2\ v (e F ) * 2 + 2x (e F ) |* 2 | 2 *2 - K (e F ) 7 (e F ) (*i - f 2 ) = 



2a 2 (e F )\^ 2 \ 2 +1; 2 (ei,)|*i| 2 - 1 = 0. 



(40) 



Solving the first two equations of [40] yields the gap equations in this regime. The solutions are 



* 2 



7 — a 

K 

x 



1 



2\ v (7 - a) 



1/2 



jo 



(41) 



where a necessary condition is that A„ be negative. Substituting these results into the final equation 
of[40]yields the transcendental equation for e F in this case, 

e {£F)l 2 (e F ) 



2a 2 (e F )\ v (ep) 



1 



K(e F )"/(e F ) 



0. (42) 



X(ep) r ' 2a 2 (e F )[ 1 (e F )-a(e F )} 2 IV ' 2A„ (e F ) [7 (e F ) - a (e F )] 
Equation [42] is solved numerically using the same algorithms and tolerances as equation 1381 for 
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Figure 2: (Color online) Plots of the Fermion/Boson number and chemical potential (ij? asa function 
of the Feshbach detuning v in the deep-BEC and shallow-BEC regimes. Here, h — m = V = 1. The 
renormalization cutoff is set to 18. The curves are evaluated for two values of the atom molecule 
coupling g r . Figures (a), (b) and (c) contain plots for the ratio iif/v, the condensate fractions and 
[if respectively for g r = 25, up = —0.3, ub = and M = 100. Figures (d), (e) and (f) contain plots 
for the ratio \ifIv, the condensate fractions and [If respectively for g r = 40, uf = —0.3, ub = 
and N = 100. In figures (b) and (e), the Fermion condensate fraction rip is shown in red, and the 
Boson molecular condensate fraction ns is shown in blue. Note that the pure molecular BEC is 
achieved when ub = 0.5, half the total number M. Also, note the onset of the shallow BEC regime 
at v ss —50 for (b) and at v ss —150 for (e). Also, note from figures (a) and (d) that fip approaches v 
for sufficiently large values of —v, indicating the onset of a pure molecular BEC in accordance with 
the results in Equation 1341 Finally, note that figs (c) and (f) show that \xf is approximately egual to 
v for large \v\, but deviates from this at smaller \v\, as can be seen more clearly in figs (a) and (d). 
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Figure 3: (Color online) Plots of the Fermion/Boson number and chemical potential fip as a function 
of the Feshbach detuning v in the deep-BEC and shallow-BEC regimes. Here, H = m = V = 1, 
up — —0.3, ub = 0.6 and Af = 100. The renormalization cutoff is set to 14.5. The curves are 
evaluated for g r — 25.0. Figure (a) contains plots for the condensate fractions as a function of the 
Feshbach detuning v. The Fermion condensate fraction np is shown in red, and the Boson molecular 
condensate fraction tib is shown in blue. Note that the pure molecular BEC is achieved when tib = 
0.5, half the total number Af. Also, note the onset of the shallow-BEC regime at v ~ —100. Figure 
(b) contains plots of the chemical potential (i^asa function of v. Note here that fj,p approaches v 
for sufficiently large values of —v, indicating the onset of a pure molecular BEC in accordance with 
the results in Equation |34[ 
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Figure 4: (Color online) Plots of Qt as a function of chemical potential (ip for g r 
ub = 0, and A/" = 100. Compare order of magnitudes with that of O m in Table [2j 



40, u F 
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representative values of the parameters. The results are shown in figures [3] Figure [3] (a) contains 
plots of the condensate fractions np (red) and ub (blue) for Af = 100, g r = 25, \up\ = 0.3, and 
ub = 0.6. Here, Ur is chosen to be 14.5. Note the shallow BEC point at v ss —100 in Figure[3](a). In 
this case, the fermion population appears to increase much faster than in the case of noninteracting 
Bosons as v is varied adiabatically past this regime. The system should populate to a full BCS state 
shortly after unitarity. 



4 Quenched Time Evolution 

Equations 1141 provide a complete description of mean field dynamics of the Timmermans' Hamilto- 
nian in equation [2j The ^s will not, in general, represent densities when the system departs from 
equilibrium. The definition of ^ in equation [T5] indicates that while |^ 2 | 2 (£) will be proportional to 
the Boson density, |\&i| 2 (t) will no longer represent the Fermion density out of equilibrium (which 
will only depend on 1 — 2cr 2 1 \I/ 2 1 2 )^ since equation [26] breaks down away from equilibrium. The gen- 
eral quasiharmonic solutions of ^1,2 have been discussed by Machida and Koyama [34], and the 
global existence of weak solutions have been established recently by Chen et al ll42ll . This paper 
focuses on the dynamics in the shallow BEC regime close to the BCS-BEC crossover as the system 
is quenched from the deep-BEC regime. 

The continued use of the single mode approximation in this case can be justified by evaluating 
the maximum characteristic time scale t q for the quench below which sizable excitations of the 
closed channel Bosons and free Fermions will occur. In order to do so, the Timmermans' Hamil- 
tonian in equation [2] is written in momentum space with the closed channel excitations included, 
yielding 

H tm = ( £k ~ L lF ) a L a k<r - ^2 a kt a -ki a k'i a k't 

k,<T kk' 

+ H (K + 2V - MB) fe q fe q + U B 5 qi 6 q^q2 5 qi 5 qi+q2,qi+qi 

q qiqiqsqi 

+ 9r X! ( 6 q a P+§t a -p+§i + h - C ) ■ (43) 
kq 

Here, the first line represents the Fermions in momentum space, with the creation (annihilation) 
operators a^ pcr (a pa ) obtained by Fourier transforms of $ a (ef> a ) from equation[2j and the free Fermion 
energies are represented by e^. The second line represents the Hamiltonian for the closed channel 



12 



quasimolecular bosons 6 q (6 q ) with energies E^ + 2v — fig. Finally, the last line represents the atom- 
molecule coupling discussed in section [2] The Fermi and Bose condensates exchange momentum 
q through this coupling. The Hamiltonian above can be written in terms of the generators of the 



SU(2) ® SU(1, 1) Lie algebra B32I1 . Assuming that most of the Bosons are in the ground state given 



by q = 0, treating b as a c-number, retaining only terms up to second order in &J. (& q ) and neglect- 
ing fluctuations about the mean field for all generators, the Hamiltonian can be diagonalized in a 



manner similar to that done for a pure Boson system in Bogoliubov theory [32], with the eigenstates 
given by the generalized SU(2) (8 SU (1,1) coherent state. If the system is in a pure quantum state, 
the energy eigenvalues are given by 13211 



E sn = (2v - n B ) \bl\+Y^ ( 2SjBk + £k " Vf) 



o )E£--(E£ + 2u-ii B +2u B \b 



(44) 



Here, Ek = \J (ek — Pf) 2 + |A + g r bo\ 2 are the energies of the Bogoliubov quasiparticles from con- 
ventional BCS theory, E^ = \/E^(E^ + 2ub|6 | 2 ) are the energies of the quasiparticle excitations 
of the pure Boson system from Bogoliubov theory [32], and sin are indices that correspond to the 
ladder operators of the SU(1, 1)/SU(2) Lie algebras respectively, with s — ±| and n e N°. If the 
energies are referenced from the ground state E_i , then the excitation energies are given by 



SE in = 2 s J^ F + \A + g r b \2 + J2(2E p + nE b p ), 

6E_ in = J2 nE l (45) 

where £k=o is taken to be 0. Evaluating equation [45] in 3D with = k 2 /2m, E® = p 2 /4m and 
a momentum cutoff at p = p c = ^ to counter ultraviolet divergences shows that the sums in 
these equations go as positive powers of the particle number density (noting that V = 1 here) 
p ~ ub|6o| 2 ~ -jk, in particular, SE_i n ~ p 5 ^ 3 . Thus, a lower bound for the excitation energies is 
given by 

E c = 2yJn 2 F + \A + g r b \ 2 . (46) 

E c ~ /O 1 / 2 (see equation I26D , which is fairly macroscopic due to the presence of the effective gap 
A+g r b a . The corresponding relaxation time t c = 1 /E c defines the time scale below which excitations 
of energy greater than E c are imparted, and excitations away from the ground state (such as closed 
channel bosons and unpaired Fermions dissociated from them) are possible. In all subsequent 
analyses, it is assumed that the quench time t q < t c and both are neglected in comparison to 
the time scales of the post quench dynamics. This effectively makes the quench an instantaneous 
impulse, and an analysis of the dynamical system based only on a Boson field and Fermion pair field 
completely captures the non-equilibrium state. 

If there is no BCS state in this system for finite values of ej?, then *i = * 2 = <]>, and the dynamics 
is equivalent to one with no atom-molecule coupling g r . Thus, equations 1141 simplify to 

^ -ia* + i/3|*| 2 * = 0, 
^ + 2i\ u ^ + 2i X \y\ 2 ^ = 0. (47) 

Substituting quasiharmonic trial solutions = ^e tQt yields 

ft* -£*# + /3|*| 2 § = 0, 
ft* + + 2x|*| 2 * = 0. (48) 

Both of these equations need to hold for two unknowns *, ft. Thus, in the absence of Boson inter- 
actions (x = 0), ft = — 2A„, and the phase of the Boson condensate fraction oscillates without any 
population transfer out of the BEC in all regimes. 

The shallow-BEC regime is one where the condensate fractions np and n B are comparable (see 
figure 0. If no a-priori assumptions are made about any of the other constants, equations [18] 
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indicate that j3 can be ignored. In the case of noninteracting Bosons, a guench from a pure BEC, 
where ^1,2(0) = (l/\/2) 9r/^F, to this regime (where the dynamics is linear) produces 

9*i 



Of 



at 



■ij (*i - - za*i = 0, 



(49) 



Substituting trial solutions *i,2(i) = ^1,2 x e int , and solving the resulting characteristic eguation 
for Q, 

1 



where 



A(e F ) 
B(e F ) 



;A(e F ) 



1±4 1 



B(e F ) 



A 2 (e F ) 



2A 



«7, 



(7 — a) (2A + K7) — K7 2 



(50) 



(51) 



Here, as always the constants in Greek letters are functions of e F . It is noted that, for sufficiently 
large Hr, B(e F ) is negative, and so fl± are both real (no damping or exponential growth). The BEC 
condensate fraction, proportional to l^l 2 ^)- oscillates with freguency fl bo = |f2 + — Substituting 
the values above, 



A(e F ) 



B(e F ) 



(52) 



For sufficiently large cutoff frequency fcjj, fif, 



A 2 (e F y 

Ur and thus it is believed to represent Rabi os 



dilations of solitons that are similar to the ones rep orted in 2004 by Andreev, Gurarie, and Radzi- 
hovsky 12611 . as well as others since 2004 B27I. l30l, bill for a quantum quenched noninteracting Fermi 
gas with a Feshbach resonance. Numerical results for fi(, are plotted for g r = 40 in figure [H 

A linear stability analysis of small displacements away from equilibrium should provide some 
qualitative insight into the long term dynamics of the system even for large displacements from 
equilibrium, such as impulse quenches. Absent the interference term *i — eqns [14] resemble the 
normal form of the Hopf bifurcation l43fl . However, this does not actually occur due to the complex 
nature of the coefficients. Writing the dynamical variables in polar coordinates as *j = a^e 1 ^ and 
dropping *i — * 2 simplifies the dynamics to yield dj = and cx t, trivial trajectories that do not 
admit to bifurcations of any kind. 

Investigating the more general cases, however, do yield interesting results. Substituting polar 
coordinates into the dynamics in eans !14l and splitting up the real and imaginary parts, results in a 
nonlinear dynamical system evolving in a 4-dimensional phase space spanned by ax 2 and <f>x 2 viz. 



a 1 

a 2 



-ja 2 sin (<p2 
K7CI2 sin (</> 2 

a — j3a\ — 7 



1 ° 2 

1 COS ( 

ai 



ai 

«2 



2A„ 



(53) 



Here, only the first of two cases discussed in the previous section are considered, where the boson 
repulsion term \ has been ignored. The radial fixed points of a 3 = can be realized for nonzero a, 
if 4>2 (t) = 4>i(t) + nir. Note that, in this limit, cos (</> 2 — <j>i) = ±1 and the rhs of the last two of eqns [53 
must be equal to each other. This is satisfied by eqn 1381 once the fixed point values of aj from 



eqns[36]are taken as a\ 



a 2 = (l-r?i/) 



with r v 



1/2 



Defining io v 



= 6n, and 



applying these expressions to the fixed points above, the system evolves at the radial fixed points 



as 



* 2 (t) 



(l-r h ,)r u e-* u >' t , 
\ u — 2«7 



(54) 



14 



Small displacements away from this orbit can be provided with complex quantities S^j as follows. 



* 2 (t) = [(1 -»?,,) r„ + e"^*. 
Substituting this into the dynamics in eans ll4l vields 

(5*1 = -i (uj u + 7 - a + 2/3r 2 ) 6V 1 - i/3r 2 ^ + i-yS$> 2 + Ci , 



(55) 



m7<5*i - i [u y + 2A + K7) <$* 2 + C 2 , 



(56) 



where Cj are constants, and terms 0(|5*j| 2 ) and higher have been dropped. Separating the system 
into real and imaginary parts yields 8±j ~ J2k 3l k & x k. where 5*1,2 = &Ei,3 + 18x2,4, and j, k = 1 — 4. 
The Jacobian matrix elements that describe this dynamics are given by Ji k = dij/dxk- The matrix 
evaluates to 



Ju 



( -[7 
[7 - a] + fit* 

7 




a] - Zfirt 




"7 



-«7 





«7 


^-3 



— «7 



«7 






/ 



(57) 



at the fixed points in egn |36[ The eigenvalues fi„ are given by the characteristic equation \ J V 
Q V I\ = 0. In this case, the equation is biquadratic and yields fi* + B u Qf, — C„ = 0, where 



0, 
B v 



~- ( I -1 vP. 

Bl ' 



1/2 



(58) 



/3r 2 2 (a - 7) + 3/3r 2 -(a- 7)' - 



2 - 3) + 1) 



2 4 , fa* - 3) (n u + 1) 

« 7 + 



fa* - ir 

(a- 7 _2r 2 /3) 2 -r^ 



X K 7 , 



k 2 7 2 + 2k 2 7 3 



a — 7 - 



2^-1 
t] v - 1 



(59) 



The quantity X>„ is the discriminant of the quartic equation. The linear dynamics is given by ll43ll 



|fe(*)>~£Lic« 



where c 3 are constants, are the eigenvectors of J v that correspond 



to eigenvalues W v , and \5x(t)) is the vector given by 8x1-4. The sign of the quartic discriminant 
T> v plays a key role in determining the nature of these trajectories. The regions of interest and the 
nature of the roots Q v are tabulated in table [1] If Q, v are pure imaginary, the trajectories correspond 
to orbital motion about the radial fixed points. The orbits constitute two oscillations of frequencies 
fij = j y/— B v /2(1 ± y/V^) 1 / 2 ]. Their combination results in periodicity iff fl^ have rational winding 
numbers n + /n~ = fi+/f2~, where are integers with no common factors. In that case, the period 
of the resultant oscillations are T v = n(n + /Q+ + n~/fl~). The changes in the number density of 



the BEC, denoted by \8^ 2 



8x1 



8x1, will have 4 oscillations viz 2f2^ and f2+ ± fl„ . The net 



period is much smaller than that of the *jS. When the Q. v start to pick up real parts, then the 
trajectories corresponding to negative real parts are stable and decay in spirals to the orbits of 
the radial fixed points, thus creating a stable limit cycle. The eigenvalues with positive real parts 
depart further away from equilibrium in an unstable limit cycle. However, these trajectories are 
not expected to actually occur in these systems since they correspond to displacements that violate 
number conservation. In general, however, the radial fixed points have a saddle-like nature in this 
regime. In the case when Q v are purely real, then the radial fixed points remain saddles, but the 
decay to the radial orbit is non-spiral and very fast. 

Note that the frequency 0, y depends on the renormalization cutoff &r through a, and so depends 
on the range of the Feshbach resonance induced interaction. For small ranges, fl u ~ kn and so this 
differs from the oscillations about equilibrium seen in BEC systems with no Feshbach resonance. 
In general, a profile of the discriminant T> v and the coefficient B v as functions of the detuning v, 
obtained using eqns [59] and [58] will give us the bifurcation diagram of this system for any choice 
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< o 


< V v < 1 


V u > 1 


B v <0 


All roots are complex 
in pairs of conjugates 
with opposing real parts 


n„ = ±y^i (i± VP.l) 

All roots are real with 
opposing signs 


/ i K I / , \ 1/2 

One pair of real roots with 
opposing signs and one pair of 
imaginary roots with opposing signs 


B v >0 


/ |H 1 / , \ I/ 2 

All roots are complex 
in pairs of conjugates 
with opposing real parts 


fl v =±i<J^ (l± y/\V v \J 

All roots are imaginary with 
opposing signs 


a, = ±iy%i(i + Vl^lj 
^ = ±^(^1-1) 1 

One pair of real roots with 
opposing signs and one pair of 
imaginary roots with opposing signs 



Table 1: Table of the eigenvalues fi„ of the Jacobian J v . The table breaks the parameter space into 
2 regions (depicted as rows) demarcated by the sign of B v . In each zone, T> v further demarcates 
to 3 sub-zones (depicted as columns). The eigenvalues, together with their natures, constitute the 
elements of the table. 

of parameters. Figure [5] (a) and (b) show plots of B v and T> v as functions of detuning for the choice 
of parameters shown in the caption. Here, B v is positive for small detuning and drops below zero 
at a critical value as \v\ is increased. Before that happens, however, V v starts from a positive value 
that is less than unity. As table Q] indicates, this means that the eigenvalues are pure imaginary. 

While B v is still positive, the discriminant drops to negative values at v « 88.6(0), and the roots 

are now complex conjugate pairs. The value of B v drops down to zero at around v = —943.77(0), 
with the roots now vanishing at the origin. As B v becomes negative, V v is also negative and the 
roots retain their ongoing structure until V v switches back to positive at around v ss —6715.(0), 
after which the roots are real with opposite signs as per table [TJ All these trends can be seen in the 
numerical evaluation of the eigenvalues plotted in the bifurcation diagram in figure [5] (c). Thus, a 
bifurcation is seen at the locus of T> v = in the shallow-BEC regime. Here the limit cycles coalesce 
into orbital trajectories around the fixed points as the detuning v is varied. However, this is not 
a Hopf bifurcation. The real parts of the eigenvalues do not cross over from the negative to the 
positive side as would be the case with a Hopf bifurcation [ 14311 . 

The paragraph above demonstrates that the 'deep-BEC ' and 'shallow-BEC regimes can be dif- 
ferentiated by the sign of the discriminant V v inside these regimes. Finally, a rapid quench from the 
former to the latter can be expected to yield orbital phase space dynamics of around the equi- 
librium fixed points after a sufficiently long time. Variations in the BEC number density, given by 
\^ 2 \ 2 — %l + x\, is thus expected to show rapid nontrivial oscillations that can qualitatively resemble 
some sort of 'collapse and revival' of the matter wave packet. 

Now, the dynamics is analyzed numerically for larger values of ep in the shallow-BEC regime. In 
this region, numerical methods are utilized to solve eauations ll4l The system is evolved numerically 
using Runge-Kutta-Fehlberg methods using the initial conditions ^2 = g^/2eF, corresponding to 
|6 | 2 = A/72 with Af taken to be 100. The quenching is assumed to have taken place by a sharp 
variation of £f (via the Feshbach detuning v) from a large negative value (in the deep-BEC regime) 
to the shallow-BEC region. Figures [6ja)-(h) show the time variation of the BEC condensate fraction 
for the two representative parameter sets from figure [2] in section [3] Figures and (e) show 
the time evolution for very small times (~ 10~ 2 units). These signals are the master signals of the 
dynamics. The master signals appear to be Rabi oscillations about the equilibrium value |# 2 | 2 with 
a frequency of Q m . In figures Ob) and (f), discrete time samples of the actual signals are plotted 
in steps of 2n/fl m . The presence of another oscillation(s) of frequency fl e indicates that the actual 
signals contain a superposition of several comparably fast Rabi oscillations interfering with each 
other, and that fl e is their beat frequency. For larger times, note from figuresOc) and (g) that these 
beats damp out over a long period of time ta, until a sudden onset of partial collapse and revival of 
the matter wave begins after a short initial relaxation time t r (indicated in the respective figures). 
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The revival time tji is also indicated in the respective figures. Approximate values of td, t r , tR and 
Sl m ,e have been obtained from figures [6] and can be seen in table [2j The numerics plotted in fig [6] 
indicate that the onset of this collapse and revival takes place at around 10 3 numerical units. The 
values of 6f , which are of the order of 10 2 in numerical units, can be obtained from table [2j Note 
that these times are expressed in numerical units. In order to obtain them in actual units, it is 
necessary to start with the system trap size L u , and note that the Fermion mass (M in actual units, 
m in numerical units) is taken to be unity in these numerical units. Thus, the relationship between 
numerical time t and actual time T can be obtained using egns 1161 1181 the dimensionless scaling 
of time by (epd)^ 1 in eguations [14] , and the fact that the base unit of time is T u = 2ML 2 u jh. The 
relationship is 



where the chemical potential ep is expressed in numerical units. Thus, the trap size L u can be ad- 
justed for a particular t until an experimentally accessible T is attained. As an example, if the atoms 
involved are Fermionic isotopes of Lithium (with mass M ~ 10~ 26 kg), and the onset of collapse and 
revival (at around t ~ 10 3 ) is desired to be of the order of milliseconds, then the equation above 
necessitates that the trap size be of the order of 300 nanometers, with the larger confinements re- 
quiring longer times T growing as VT. Ultracold BCS systems with such dimensions have already 
been obtained using superconducting atom chips [44], where trap lifetimes can range from seconds 
to minutes B45I1 . Therefore, such time scales are experimentally accessible. The similarities in order 
of magnitude with fl m from tableland Q& from figure[4]suggest that the collapse and revival effect 
may be soliton-like Rabi oscillations reported in the literature. Also, note from these figures that 
all three time scales, td, t r and tR, are larger for smaller g r . Thus, it is expected that td — >• oo as 
g r — > 0, removing the collapse and revival effect. The master signal frequency is also expected to 
reduce to 2A„ in this limit. It is also noted from figureEtd) that the revival process characterized by 
tR seems to be interrupted by a signal that has the same temporal characteristics as the relaxation 
signal characterized by t r , thereby inhibiting a total revival of the matter wave. 

In the case of interacting Bosons, the parameters from figure [3] are used and the coefficient /3 
neglected to simplify the dynamics. Thus, equations 1141 simplify to 



This system is evolved numerically using the same initial conditions and numerical methods as 
detailed in the previous paragraph (however, the tolerance required for convergence was an order 
of magnitude lower). The Boson interaction amplitude was taken to be fairly small so as to see 
the effects of departing slightly from the ub = case for equations 1491 For small values of ub, 
only a damped oscillation was seen for all the time intervals that were allowed by the numerical 
tolerances. For ub = 0.6, the collapse and revival effect was seen for large times, just as it was with 
the case of noninteracting bosons. 

In both of these cases, the nonlinear contribution to the dynamics go as \x\ 2 x, absent which the 
dynamics is oscillatory. Therefore, it is believed that this effect is caused by the nonlinear dynamics 
from this type of expression alone. The results of the numerical simulation for ub = 0.6 are shown 
in figures [3a) through (d). The frequencies Q m and fi e , as well as the time scales td, t r and tR, are 
tabulated in table [2j Note that the decay time td is much higher for the case of interacting bosons 
and no /3, suggesting that the contribution of the Boson interaction to the collapse and revival is 
much smaller than that of the term that is nonlinear in the gap parameter and weighed by /?. These 
striking phenomena should manifest in the time-of- flight absorption images of the gas. 

5 Conclusion 

In this paper, the dynamics of an ultracold gas of Fermions with a narrow Feshbach resonance 
has been formulated in the dual-channel case after making the single mode approximation for 
Bosons. The mean field nonlinear and complex TDGL dynamics of the Fermions, coupled to the 
Gross Pitaevski dynamics for the composite Bosons, have been obtained. The Fermion dynamics is 




(60) 



*i + i7 (*i - \& 2 ) - ict^i = 
*2 + 2a* 2 + 2ix|*2| 2 *2-iK7(*i-*2) = 0, 



(61) 
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g r -v e F u B 




n e 


t d 




tR 


25 55 125.78 


1400 


10 


400 


9 


6.0 


40 140 192.40 


3500 


6 


125 


4 


2.5 


25 110 85.5 0.6 


1050 


1 


600 


12 


8.0 



Table 2: Table of master and envelope frequencies Vl m>e , as well as decay, relaxation and revival 
times (td, r ,R respectively) of the quenched dynamics shown in fiqures[6]and[7l The system is evolved 
for two values of g r . The correspondinq vs and e^s are also shown. The values of the parameters 
are chosen as a representative sample, based on our choice of lenqth and atom number scales as 
explained in section [3j and the conclusions are not affected qualitatively if the ranqes are varied 
within these orders of maqnitude and kept in the reqions of interest. 



encapsulated in the temporal variation of the superfluid qap parameter A, and that of the Bosons 
in the qround state amplitude b . 

The equilibrium states for the system in the BEC side have been evaluated by lookinq at the 
adiabatic evolution of the stationary solutions of the dynamics, which leads to a population transfer 
from the BEC superfluid to the BCS superfluid for nondiverqent Fermion interactions. The dynamics 
of the system as it is quenched from a pure BEC state to a state in the shallow-BEC reqime (accom- 
plished by a rapid variation in the Feshbach detuninq v) has been analyzed numerically. Nonlinear- 
ities in the dynamics cause the Rabi oscillations to relax back to it's equilibrium state. However, 
at larqe times, the relaxation qets interrupted and a collapse and revival type phenomenon of the 
Bosonic matter wave field ensues. This effect seems to be caused by the interference between the 
multiple modes of the dynamics, and shows a strikinq phenomenon that is analoqous to the collapse 
and revival effect that has been reported experimentally for Bosons. 
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Figure 5: (Color online) Bifurcation diagram in the complex plane of Q v obtained by numerically 
solving the BCS gap and number equations, applying the results to the Jacobian, and diagonalizing 
it. The numerics are performed for 100 particles with g r = 25, up = — 0-3, ub = 0, and renormalization 
cutoff k r = 18.0 in a unit volume for unit particle mass, where h has also been set to unity. Figure 

(a) shows the coefficient B defined in eq [59] as a function of detuning v. Note that B > until 
v ps —943.77(0) and then switches sign. Figure (b) shows the magnitude of the discriminant T> as a 
function of detuning v in a semi-log plot. The regions of positive (negative) V are indicated in blue 
(red). Figure (c) shows the actual bifurcation diagram, with the roots Q shown in the complex plane. 
The origin, where B = (see table [TJ, is indicated, as well as the points where the discriminant V 
vanishes. The latter occurs at v ~ —88.6(0) and v ~ —6715.(0). The arrows indicate the direction 
where the magnitude of the detuning v (given by \v\) increases in the range shown in figs (a) and 

(b) . 
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Figure 6: (Color online) Plots of the signal | \I>2 1 2 ) from the solutions to eguations [14] for ug = 
after the system is guenched from a pure BEC to the shallow-BEC regime (see figure [2). Here, 
h = m = V = 1, and up = —0.3. Figures (a) - (d) show plots of the time evolution for g r = 25.0 
and v = —55.0. Figures (e) - (h) show plots for g r = 40.0 and v = —140.0. Figures (a) and (e) show 
plots for very small times (~ 10~ 2 ). Figures (b),(c) and (f),(g) show plots of discrete time samples of 
the signal taken in units of the local time period from figures (a) and (e) respectively. Thus, figures 
(b),(c) and (f),(g) show the signals that envelope the master signal in (a) and (e) respectively For 
figure (a), the master signal has frequency Cl m ~ 1400, and for figure (e) Cl m ~ 3500. Figures (b),(c) 
and (f),(g) show the same corresponding envelope plots for two different time intervals. Partial 
collapse and revival of the matter wave can be seen after a lengthy decay time, where the original 
envelope damps out as it oscillates. The envelope oscillation freguency fl e ~ 10 for figure (b),(c) 
and ~ 6 for figures (f),(g). The decay times td for figures (c) and (g) are approximately 400 and 125 
respectively. The region where collapse and revival take place is magnified and shown in figures 
(d) and (h), where the entire signal (not just the envelope) is now being plotted. The presence of 
partial collapse and revival of the matter wave can be noted after a relaxation time t r as indicated 
in figures (d) and (h). The revival times i# are also indicated in these figures. In figure (d), t r ~ 9, 
t R ~ 6. In figure (h), t r ~ 4, t R ~ 2.5 
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Figure 7: (Color online) Plots of the signal |\&2| 2 (i) from the solutions to eguations [49] for ub = 0.6 
after the system is guenched from a pure BEC to the shallow-BEC regime (see figure [2}. Here, 
h = m = V = 1, Uf = —0.3, and g r = 25. Figures (a) shows plots for very small times (~ 10~ 2 ). 
Figure (b) shows plots of discrete time samples of the signal taken in units of the local time period 
from figure (a) . Thus, figure (b) shows the signal that envelopes the master signal in (a). For figure 
(a), the master signal has freguency fl m ~ 1050. Figures (b) and (c) show the same corresponding 
envelope plots for two different time intervals. Partial collapse and revival of the matter wave 
can be seen after a lengthy decay time, where the original envelope damps out as it oscillates. 
The envelope oscillation freguency fl e ~ 1. The region where collapse and revival take place is 
magnified and shown in figure (d). The presence of partial collapse and revival of the matter wave 
can be noted after a relaxation time t r as indicated in figure (d). The revival times are also 
indicated in these figures. In figure (d), t r ~ 12, t R ~ 8. 
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